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The nucleon self-energy and its relation to the nucleon axial charge g A are discussed 
<^ at large N c . The energy is compared for the hedgehog, conventional, and recently 

proposed dichotomous nucleon wavefunctions which give different values for gA- 
We consider their energies at both perturbative and non-perturbative levels. In 
perturbative estimates, we take into account the pion exchanges among quarks up 
43 to the third orders of axial charge vertices, including the many-body forces such as 

the Wess-Zumino terms. It turns out that the perturbative pion exchanges among 
Qh valence quarks give the same leading iV c contributions for three wavefunctions, while 

their mass differences are O(Aqcd)- The signs of splittings flip for different orders of 
the axial charge vertices, so it is hard to conclude which one is the most energetically 
favored. For non-perturbative estimates involving the modification of quark bases, 
J> we use the chiral quark soliton model as an illustration. With the hedgehog quark 

wavefunctions with gA of 0(N C ), we investigate whether solutions with coherent 
pions are energetically favored. Again it is hard to give decisive conclusions, but it 
is likely that adding the confining effects disfavors the solution with the coherent 
pions, making a pion cloud around a nucleon quantum rather than coherent. The 
nuclear matter at large N c is also discussed in light of the value of gA- 
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1 Introduction 



Recently we discussed the phenomenological problems for nucleons at large 
N c , which are caused by the coherent pions surrounding nucleons [1]. The 
source of a coherent pion cloud is the large nucleon axial charge gA of 0(N C ), 
whose value depends upon quark wavefunctions inside a nucleon. 

The problems related to large gA or coherent pions manifestly appear in the 
pion exchange part of the nuclear potential. Consider dilute nuclear matter 
where the pion exchange is the most important. For the nucleon-nucleon po- 
tential, the overall strength of the pion exchange is proportional to the square 
of the vertex, (g A / f w ) 2 ~ N c . The kinetic energy of nucleons are at most 
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~ jP/Mn ~ Aqcd/^c, so the kinetic energy is much smaller than the po- 
tential energy. It means that at large N c , nucleons localize at the potential 
minima, forming a crystal even at dilute regime [2,3,4,5]. In reality, however, 
dilute nuclear matter is like a liquid where the potential and kinematic effects 
are equally important, and its binding energy is 0(1 — 10) MeV, hardly recog- 
nizable from the QCD scale Aqcd ~ 200 MeV. For this qualitative discrepancy 
from the N c = 3 world, usually the large N c results are not regarded as useful 
guidelines for studies of a nuclear matter. 

This situation, however, might be altered if the value of qa were not 0(N C ) but 
0(1). With qa ~ 1, both the potential and kinetic energy become comparable, 
~ Aqcd/N c ~ 60 MeV at N c = 3, so it is much easier to expect the 0(1 — 10) 
MeV binding energy, by cancelling out the energetic cost and gain at the same 
order. Then the large binding energy problem could be solved, or at least be 
largely reduced. In this picture of small qa , the long distance part of nucleon 
potentials should be described not by coherent pions but by a few quantum 
pions, and in fact, the latter has been a description in conventional nuclear 
physics [6,7,8]. Motivated by this observation, in Ref.[l] we started our attempt 
to construct the ground state nucleon wavefunction with qa ~ 1 in case of the 
SU(N C ) gauge theory. In this paper, we shall continue our discussions about 
a nucleon wavefunction with small qa- 

This reduction of qa itself does not solve all of the problems, however. We 
must handle a large attractive force of 0(iV c ) at intermediate distance as well. 
It comes from the exchange of the a-meson whose mass is supposed to be 
sensitive to N c in contrast to other mesons [9,10]. We imagine that at large 
N c , the range of o is comparable to the range of the strong repulsive force 
originating from the u meson, and the spatial size of the attractive pocket 
becomes very small. Then the binding energies of nucleons would arise not 
from the leading N c contributions, but from the l/N c corrections. This aspect 
will be discussed elsewhere, but in what follows, the longest range part of the 
interaction must be first resolved before arguing the forces at intermediate 
distance. Therefore, as a first step, we focus on the pion problem. 

The source of the axial charge qa is quark dynamics inside a nucleon. In this 
work, we will consider the valence quark contributions to the axial charge, 
g A al , that is estimated from the valence nucleon wavefunction. As a theoretical 
guideline, in the following we consider problems in the context of the con- 
stituent quark picture, including the meson exchanges. If we include only the 
pions, the description we shall use is analogous to the frameworks of Manohar 
and Georgi [11], Weinberg [12], and Glozman and Riska [13]. Following a study 
of Weinberg [14], the quark axial charge is assumed to be 1. 

For considerations of the axial charge, it is convenient to use the SU (4) spin- 
flavor generators [15,16], formed by isospin, spin and axial charge operators: 
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(i) 




where g means the operator acts on the g-th quark. We will not require that 
the Hamiltonian is SU (4) symmetric, and in fact, nucleons are not eigenstates 
of axial charge generators. Nevertheless the use of the £77(4) makes it easier 
to get several qualitative insights. For instance, the diagonal axial charge of 
the proton with spin up is given by 



for the valence quark part. When we write \pt), h means the valence quark 
part for which we do not include extra contributions from induced pion clouds, 
etc., unless otherwise stated. 

Below we will consider three types of the valence baryon wavefunctions in 
order. The hedgehog, conventional, and recently proposed dichotomous wave- 
functions [1]. They provide different values of g™ x : —N c , ±(iV c + 2)/3, and ±1, 
respectively. 

In many studies, large N c baryons are discussed starting with the baryon wave- 
function of the hedgehog type [17]. The hedgehog state is a mixture of different 
spins and isospins but has a definite and very large axial charge, g™ 1 = — N c , 
that generates a large pion field. It allows us to apply the mean field or coher- 
ent picture of pions which is an underlying basis for the Skyrme type models 
[18,19,20], its holographic version [21,22,23], or chiral quark soliton models 
[24,25,26,27,28]. The states with definite isospin and spin can be obtained by 
collectively rotating the hedgehog wavefunction and the pion cloud. Therefore 
nucleons constructed in this way (we call them hedgehog nucleons) accompany 
a coherent pion cloud. 

A more theoretically sound way to compute nucleon properties is to always 
keep the quantum numbers during calculations, although such computations 
are technically much more involved. It is nontrivial whether such a nucleon 
involves a coherent pion cloud or not. This is because once we fix the isospin 
and spin of a nucleon state, it is inevitably a mixture of different axial charge 
states (See discussions later). 

In the conventional nucleon wavefunction, all quarks occupy the lowest £- 
wave orbit, and the spin-flavor wavefunction is maximally symmetrized to 
satisfy the Fermi statistics. The conventional nucleon wavefunction produces 
g™ 1 = ±(iV c + 2)/3 at odd N c , while at even N c it becomes zero [1] (here the 
"nucleon" at even N c simply means the lowest isospin and spin state). This big 
sensitivity of g™ 1 to whether N c is odd or even is a consequence that a nucleon 
wavefunction is a superposition of different axial charge states. In particular, 



9a = <Pt|i2ssbt>, 



(2) 
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for an even N c "nucleon" , the contributions from different axial charge states 
completely cancel out one another. 

In the dichotomous wavefunction at odd N c , we place N c — 1 quarks in the 
lowest ,S-wave orbit in such a way that they have zero isospin and spin. Then 
the remaining quark carries the same quantum number as nucleons, and we 
place it in a spatial orbit different from the other N c — 1 quarks. As seen from 
the computation for the N c even "nucleon", the contributions to g™ 1 from 
the N c — 1 quarks cancel out one another, and g™ 1 is saturated solely by a 
single quark. It means that g™ 1 is 1. An obvious question, however, is how 
to assure that the dichotomous state is energetically favored compared to the 
conventional one. To place one quark in an orbit different from the lowest 
S- wave orbit, we have to pay an energy penalty of O(Aqcd)- 

Obviously which states appear as the ground state depends on the quark 
dynamics inside of the baryons. The purpose of this paper is to give the argu- 
ments on the energetic differences. Unfortunately, it turns out that it is very 
hard to derive a definite conclusion about which one is the lightest state. 

We will argue the energy dependence of a nucleon on in perturbative and 
non-perturbative contexts. 

In a perturbative context, we consider pion exchanges among quarks inside 
a nucleon, which are supposed to be strongly correlated to the axial charge 
of a nucleon. The summation of their contributions depend on the spin-flavor 
wavefunctions, and are computed group theoretically. We evaluated the self- 
energy of baryons up to the terms quadratic and cubic in the axial charge 
operators. The Wess-Zumino vertices [29,30,31] are included. 

We will argue that at the level of peturbative evaluations, the leading N c 
contributions for hedgehog, conventional, and dichotomous wavefunctions are 
energetically degenerate. The pion exchange contributions for these wavefunc- 
tions differ by at most 0(Aq CD ), at each order of the axial charge operators. 
So three wavefunctions differ in energy by O(Aqcd) unless there are subtle 
cancellations among different orders of axial charge operators. 

The computations based upon the perturbative pions are certainly insufficient, 
especially when we argue the nucleon with a large axial charge. When the 
strength of the pion field becomes very large, such a large pion field may drill 
a hole in a media of the chiral condensate, costing energy. Accordingly the 
constituent quark wavefunctions are modified as well. In this regime, we must 
take the non-perturbative effects into account. 

This kind of contribution is considered in the chiral quark soliton model in- 
cluding a full evaluation of the quark determinant [25,26,27,28], provided that 
the background coherent pion field is of the hedgehog type. In this picture, 
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the production of the coherent pion cloud costs energy of 0(N C ), while quarks 
acquire energetic gains of 0(N C ), by being bound to the pion cloud. The opti- 
mized configuration of pions and quarks is determined by the energetic balance 
between these 0(N C ) contributions. 

We will see that for the hedgehog wavefunctions of quarks, it is justified to 
take the pion fields to be static and classical. Within a context of the chiral 
quark soliton picture, we give arguments why we expect that the hedgehog 
wavefunction is not energetically favored. On the other hand, for the conven- 
tional or dichotomous nucleon wavefunctions, we do not know how to take the 
appropriate mean field for pions, so we could not show the nonperturbative 
expressions of energies for these two wavefunctions. 

In this paper, we discuss only two flavor cases, close to the chiral limit. Our 
metric is g^ u = diag(l, — 1, — 1, —1). 

This paper is organized as follows. Sec. 2 is devoted to some preparation of 
basics which are used in later sections. We start with counting of the quark- 
meson coupling and classify its total contribution based on the spin-flavor 
wavefunctions of nucleons. Then we briefly review basic relations of the SU (4) 
generators to the extent necessary for our arguments. In Sec. 3, we show the 
construction of the hedgehog, conventional, and dichotomous wavefunctions, 
and compute their axial charges. In Sec. 4, we compare the energies of three 
wavefunctions at the perturbative level. In Sec. 5, we give non-perturbative 
considerations within a context of the chiral quark soliton picture. Sec. 6 is 
devoted to summary. 



2 Meson exchanges for the baryon self-energy: Preparations 

In this section, we first argue how the meson exchanges among quarks con- 
tribute to the baryon mass. Especially we will focus on the meson exchange 
in the pion channel, because its quark-meson vertex has a form of the axial 
charge operator. For evaluations of the meson exchange contributions, it is 
useful to employ the SU{A) spin-flavor algebra and its Casimir values at vari- 
ous orders. The many-body forces, including the Wess-Zumino vertex, will be 
discussed as well. 



2.1 N c - counting 

Suppose diagrams in which gluons are exchanged between two quarks. Di- 
agrams that are obtained by iterations of a single gluon exchange may be 
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a) 1 2 b ) 




Fig. 1. An example of two particle irreducible diagrams in which all gluon lines are 
crossed. We attached numbers to vertices to identify their locations after deform- 
ing quark lines, (a) Two gluon exchanges, (b) A deformed version of the diagram 
(a) which can be interpreted as a quark exchange, (c) A color line representation. 
There is a closed color line which provides an enhancement factor N c . (d) The 
corresponding meson exchange between quarks. 

treated within the one gluon exchange potential. On the other hand, there 
are also diagrams in which gluon exchange lines are all crossed, and are two 
particle irreducible (Fig.l). For the latter, by deforming quark lines it becomes 
manifest that the diagram is a planar one, and it is natural to interpret this 
quark exchange diagram as a meson exchange between two quarks. 

Here one might think that at a distance scale among quarks inside of a baryon, 
it is not appropriate to consider quark exchanges as the meson exchanges. 
Such a terminology should not be taken too literally in this paper. The word 
"meson" is used as a representative of the two particle irreducible diagram 
such as Fig.l (a), which carries the same quantum number as mesons. The 
person who dislikes this terminology can directly resum diagrams within the 
language of quarks and gluons. The real problem we can not rigorously handle 
in this paper is rather the higher towers of mesons appearing at the distance 
scale ~ Aqq D , although it is likely that such mesons do not alter the main 
lines of our thoughts. 

The meson exchanges among quarks can appear in the 0(N C ) contributions 
for the baryon self-energy, because there are many possible ways to choose 
the quark lines. We must keep them in the computation of the baryon mass 
as the leading order contributions. In the language of the effective Lagrangian 
for mesons and baryons, the meson exchanges among quarks correspond to 
the meson loops attached to a baryon line (Fig.2). 

As a reminder of the self-energy evaluation, consider diagrams such that one 
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Fig. 2. The self-energy diagram for a baryon. (a) The meson exchange among quarks. 
There are many possible ways to choose quark lines, ~ A 2 , whose contributions must 
be summed up. (b) The corresponding diagram in which a meson loop is attached 
to a baryon line. 



T * 




Fig. 3. An example of the many-body forces in the baryon, induced by mesons. The 
digaram can appear at the leading iV c contribution. 

meson is exchanged between two quarks. Following the usual A^-counting, the 
intrinsic quark-meson coupling is given by g qq u ~ N' 1 ' 2 . One can choose two 
quark lines out of N c quark lines, so its combinatorial factor is ~ A" c 2 . So if all 
of contributions are simply added without cancellations, one can estimate the 
contributions from these processes as (g qq M) 2 x A" c 2 ~ N c . 

Actually, whether the total of the meson exchange contributions are large or 
not, strongly depends upon a spin-flavor wavefunction of a baryon. 

The simplest example is the charge-charge interaction mediated by the uq- 
meson exchange (the zero-th component of the u meson). The Uq couples to a 
quark number. Since all quarks have quark number +1, its sign always appears 
in the same sign. It means that all contributions are simply additive, so the 
self-energy is {g qqM f x A" 2 ~ A" c . 

There are many other examples in which various contributions cancel one 
another. Consider the po exchanges which are responsible for the interactions 
among isospin charges. In contrast to the quark number, isospins of quarks 
may be either positive or negative, so that their contributions can cancel one 
another. For baryons with isospins of 0(1), the po -exc hange diagrams happen 
to cancel out one another, producing negligible contributions, g^ qM x 1 ~ 1/N C . 
This process is important only for baryons with isospins of 0(N C ). 



7 



In addition to the meson exchanges between two quarks, one must take into 
account many-body forces as well, because a quark density inside of a baryon 
is not dilute and there are many chances many mesons meet at once (Fig. 3). 
Although the ri-meson vertex behaves as ~ iV c 1 ~ n/ ' 2 and is suppressed, there 
are many ways to choose n-quark lines out of iV c -quark lines, ~ N™. Each 
quark- meson vertex has ~ iV" 1 / 2 , so as a total, the contributions from many- 
body forces can be ~ 

N i-n/2 x N n x N -n/2 ^ /v" c . In general, the tree dia exams 
of meson lines can appear at the leading order of N c . 

This sort of arguments are also useful for the classification of the two-body and 
many-body forces in a baryonic matter. In fact, it has been used to interpret 
the hierarchy of the several meson exchange forces in the nucleon-nucleon 
potential [32]. 



2.2 The axial charge operator in the SU(4) spin-flavor algebra 



In this work, we will mainly focus on the diagrams involving pions, because we 
expect them most directly related to the axial charge of a baryon. Considering 
the non-linear realization of pions and quarks, their couplings contain at least 
one derivative. The lowest order coupling is given b}f*~| 

Zqqir = 3 A Ql^aQ ■ (3) 

We will take the quark axial charge g\ to be 1 [14]. 

We will take the non-relativistic approximation for quarks. This approximation 
is not rigorous for the quantitative estimates for the light quark system. But 
we expect that the relativistic corrections change overall strengths of the pion 
interactions equally for different wavefunctions, without affecting their energy 
ordering. 

Applying the non-relativistic approximation for quarks in the baryons, the 
above vertex is simplified to 

(B\q W r a q(x)\B') — ► ^ J2{B\r^a^S(x - r<^\B') , (4) 



where the zero-th component vanishes at the non-relativistic limit, and dom- 
inant contributions come from the spatial components. This matrix element 
carries the wavefunction dependence of the pion-baryon coupling. We shall 
evalute it using the group theory. 



1 Our convention is (O^^^-^a;)!^^)) = -i^/^e iqx , U(x) = e 171 "" 7 "/^, for 
which the pion decay constant is f n ~ 93 MeV. 
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The axial charge operators, together with isospin and spin operators, form the 
SU (4) algebra. Writing operators 



r a = J2 ri q) , <n = E °P ,R« = E = E tP*P , (5) 
q i q q 

their commutation relations are 

[7~a,n] = ZitabcTc = ^ijk^k , 

[ T ai Rbj ] = ZitabcRcj , Rbj ] = ^ijkRbk , 

[ r , Oj ] = , [ R ai , R bj ] = 2ie abc 5ijT c + 2ie ijk 5 ab a k . (6) 

Clearly, the operators R ai transform as vectors under the operations of isospin 
or spin rotations. 

In the following, we classify the baryon states by specifying their eigenvalues 
of the Cartan bases for the SU (4) algebra. As Cartan bases, we choose 

( r 3 , a 3 , #33 ) , (7) 

where these operators commute one another so that we can label states with 
these quanta. The states in the irreducible representation V can be labeled as 



Irg,^,^;^). (8) 

For instance, fundamental representations are formed by single quarks 

H> = |i, l, l) , |dt> = |-i, i, -i) , 1*4) = |i, -i, -i) , \dl) = |-i, -i, i) , 

(9) 

which will be building blocks for the baryon wavef unctions. We will give their 
explicit forms in the next section. 

2.3 The processes involving pions and their relations to the axial charge 



Now let us see how these operators appear in the diagrams involving pions. 
We classify contributions by the order of the R operators. For instance, one 
pion exchange between two quarks is the order of R 2 , the diagram involving 
four pion emissions from four quarks and one 4-7T vertex is the order i? 4 , etc. 
Later we will evaluate the expectation values of these generators. 
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Fig. 4. The one pion exchange contribution to the baryon mass. Its total size is 
proportional to the baryon matrix element for the R 2 ai operator. 

2.3.1 The R 2 operators: One pion exchange 

The simplest diagram contributing to the R 2 operator is the one pion exchange 
among quarks. We wish to evaluate the following matrix element (Fig.4) 

i 2 f 

M 2 n = ( B \^pJ dxd V (diTTa(x) qia*,T a q{x)) (dj-Kbiy) qijl5nq(y)) \B) , (10) 

where we dropped off the zero-th component of the axial- vector current, taking 
the non-relativistic limit, Eq.(4). M.2-K should be matched with the potential 
as M 27T = -i(Vto)T, because ( e ~'< Ho+v ^ T ). 

We shall carry out the integration of virtual lines, leaving only quark operators. 
The pion propagator is 

(TdMx) dMv)) = s ab j^ e-*-»> i . (11) 

In the non-relativistic limit, p\ term in the denominator is smaller than the 
spatial component by \p\/M where M is the consistuent quark mas 50 Taking 
the non-relativistic limit, we decompose the propagator into the central and 
tensor parts, 

PiPj ) 5 a f flW p 2 

p2 _ m 2 3 ^2 _|_ m 2 y p 2 3 J p 2 + ml 

The second term is activated in the intermediate processes such that two 
quarks in the S- wave are scattered into the D-wave. We shall ignore such a 
contribution, and will take into account only the first term. Together with 



2 This can be seen by writing equations in terms of the old-fashoned perturbation 
theory. 
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these simplifications, 



(Td<K a (x)dMv)) -^^o-»)) ^-i?)-^ ]■ (13) 

The reason why the 5-function appears is that the derivative coupling of pions 
to quarks becomes larger for larger momentum transfer. It cancels out the 1/p 2 
behavior of the pion propagator, so that the pion propagations with relatively 
high momenta are not suppressed (although in a more realistic treatment there 
should be the UV cutoff due to form factors). On the other hand, the second 
term is proportional to the pion mass, so disappears in the chiral limit. 

In the chiral limit together with aforementioned approximations, Eq.(lO) is 
reduced to 

Eq.(lO) ->■ J d 4 x (B\qw 5 T a q(x) q-fa 5 T a q(x)\B) 
= 7^/d 3 :?EE< 5 l^^ (14) 

iZ J r 7r J q , q ' B' 

We mainly consider all quarks in \B) occupy the same spatial orbit, ip. To 
have nonzero matrix elements, quarks in \B') must occupy ip except the g-th 
quark which hits the (^-function. We ignore \B') including a radial excitation 
of a quark, and consider only \B') with all quarks in the spatial orbit ip. Then 
the equation is further simplified by summing the intermediate state \B'), 

Eq.(14) ~ ^ /d 3 f |y,(f)| 4 (B\Rl\B) SF . (15) 

where \B) SF means the spin-flavor part of \B). Therefore the potential is 

<v- 27r > ^ - y i^ ( ^)i 4 <^i^i^> sf — 4It^ <^i^li^>s F , (i6) 

which reduces the baryon mass. Here we assumed that the spatial extent of 
ip(x) is ~ Aqq D , so that |y?(:?)| 2 ~ Aq CD from the normalization. 



2.3.2 The R? operators: Wess-Zumino vertex 

Similar manupulations can be done also for higher order products of the R 
operators. Here of particular interest is the u — 3ir vertex coming from the 
Wess-Zumino term (Fig. 5). The diagrams involving the u meson are very 
important because u strongly couples to a baryon due to its large quark 
number density of 0(N C ). In fact, in the Skyrme or chiral soliton models, this 
term is responsible for a repulsive force of 0(N C ) between a quark number 
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Fig. 5. The oj — 37r contribution to the baryon mass, mediated by the Wess-Zu- 
mino vertex. Its total size is proportional to the baryon matrix element for the 
e a bc£ijkRaiRbjRck operator. 

density and a topological density, which stabilizes the configuration of the 
chiral soliton. 

The quark-o; meson coupling appears in the form 

£u = -g u Un qi^q ■ (17) 

If we take the sign of g u in this way, the Wess-Zumino term for the u — 3tt 
coupling is given byp~| 

= T^f^ e^e abc w„ (cV a ) (d a n b ) (d fi n c ) . (19) 

When evaluating this contribution for the baryon states, the term involving 
ido field is much larger than terms with aVs. Thus we consider only 

£a)0-37T = - JT^fk e ijk£abc w (cV a ) (<9j7T 6 ) (d k n c ) . (20) 

As before, we draw a diagram and integrate the virtual meson propagations, 
leaving the matrix element for quark operators. We compute 

3 



J d 4 xd 4 yd 4 zd 4 wd\(B\r w _ 3n \B) , (21) 



_i N ^ 



3 The expression can be intuitively checked from the cj-quark vector coupling, 
£qqu = — ff^w^J^, with the Goldstone-Wilczek's topological quark number current, 

JM = "2^2 ^ Q/3 t r/ [(rfd v U)(Utd a U)(Utdf,uj\ , (18) 
where the trace is taken over the SU(2) fundamental representations. 
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where 



1 w — 3-rr (g7i75r a g(a:)) (qijl$nq(y)) (qikl^ c q(z)) (qioq{w)) 
x (<9;7r a (:r)) (dj-K^y)) (d k 7i c (zfj u (w) 

x e /mn e de/ u (u) (dtn d (u)^ (<9 m 7r e (w)) {d n ir f (u)^ . (22) 
The propagator of the u meson, in the non-relativistic approximation, is 

-mn 5(wq - u ) G w (w - u) . (23) 



where 



p— m u | to— u\ 

G^{w-u)= (24) 

47T IU — U\ 



Assembling these approximations together, our vertex now takes the form 
N a 2 T r 

M u -to^i ^Wf* J d x d ™( B \ R{3) (x) G UZ-™)Qo(w)\B), (25) 
where we defined the cubic product of the axial charge operator, 

RW(x) = e abc e ijk R ai R bj R ck (x), R ai (x) = £ 5(x - r<*>) , (26) 



and the quark number operator is Qq{w) = qioq(w). The quark number is 
very large, N c , for any baryon states. Again assuming that all quarks occupy 
the spatial orbit </?, then Eq.(25) becomes 



NaT r 
^0^r (B\R^\B) SF J dxdw \<p{x)\*G u {x-w) x (NMw)f) . (27) 



This takes the form of the interaction between the distribution of the R 3 
operator and the quark number density. Finally the potential energy is 

= ~ ^33^6 ( B \R^\B) SF J dxdw \<p(x)\°G„(x-w) \ V (w)\ 2 . (28) 

As we will see, the R 3 term appears to be negative and 0(N 3 ). The overall 
factor is N 2 g 2 /f® ~ l/N 2 , so K>-3tt gives the energetic cost of 0(N C ). 
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3 Hedgehog, conventional, and dichotomous nucleon wavefunctions 
and their axial charges 



In this section, we show explicit forms of the hedgehog, conventional, and 
dichotomous wavefunctions. Then we compute the expectation value of the 
i?33 operator. 



3.1 A hedgehog nucleon wavefunction 



The baryon wavefunction giving a hedgehog type configuration has zero grand- 
spin, G = r + o = 0. The basic building block is a single quark wavefunction 



(29) 



which vanishes under the operation G = f + a. (This choice is unique. Other 
wavefunction such as (tti +d"\) does not vanish under the operation of 71 + 01, 
etc.) The hedgehog baryon wavefunction is 

im-^ff. (30) 

This is just a product of the same single quark wavefunction, so it is maximally 
symmetric in spin-flavor. The special property is that \H) has a definite R33 
eigenvalue. Indeed, in the SU(4) bases (|r 3 , 03, -R33)) one can write 

(ui -dt) = |1, -1, -1) - I - 1, 1, -1) , (31) 

where the eigenvalue -R33 is common, so that its any products have a definite 
i? 33 eigenvalue. But they are a linear superposition of the different (r 3 , 03). In 
these bases, the hedgehog wavefunction is written as 

l^> = ( ^ J E y/NcOn (-l) n |JV C - 2n, -(iV c - 2n), -iV c ) MS , (32) 

where | ■ ■ -)ms means that the state is maximally symmetrized. Here a factor 
a/ N c C n is included to normalize the state as (r 3 , 03, -R33I73, 03, -R33)ms = l PH 



For instance, 

(u| -(it) 2 = (u|) 2 + (d|) 2 - («4 dt +d| «4-) 

= |2, -2, -2) MS + 1-2,2, -2) MS - V2 |0, 0, -2) MS • (33) 
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All the states contained in the hedgehog wavefunction have large R33 eigen- 
values of 0(N C ). So fluctuations of the R33 charge, that are caused by the 
577(4) breaking terms in the Hamiltonian and happen in the time evolution, 
can be regarded as small fluctuations around the big mean value of 0(N C ). 
Thus the axial charge of the hedgehog state may be treated as classical and 
static mean field. We will use this observation when we give a nonperturbative 
consideration for the baryon self-energy. 

The axial charge of the hedgehog state can be computed readily, and is given 
by (H\R 33 \H) = — iV c , hence = —N c . The minimal eigenvalue i? 33 for N c 
quarks are —N c , so the hedgehog state saturates the lowerbound of the baryon 
axial charge g\. 

3.2 A conventional nucleon wavefunction 

The conventional nucleon wavefunctions have definite isospins and spins, while 
they are superpositions of different axial charge states. To construct states 
with small isospins and spins, it is convenient to prepare a \I — 0, S — 0) pair 
wavefunction as a building block: 



which is a superposition of pair wavefunctions with different axial charge 
states, -R33 = ±2. 

Before going futher, we stress that the construction of a nucleon from diquarks 
has nothing to do with any dynamical assumptions such as diquark correla- 
tions. In fact, at large N c the diquark correlation is suppressed. Here we simply 
use the diquarks just for the convenience to prepare appropriate spin-flavor 
quantum numbers. 

By just taking a direct product of these states, we can obtain spin-flavor func- 
tions with I — S — 0. To construct N c odd nucleons, we must further add 
one more quark to achieve the color singletness and appropriate spins. Then 
we totally symmetrize the spin-flavor-space wavefunctions for products of di- 
quarks and an unpaired quark. Usually a nucleon with a maximally symmetric 
spin-flavor wavefunction are supposed to be the lowest energy state, because 
in this way one can place all of N c quarks in the S'-wave spatial orbit relative 
to the center of the nucleon. 





(34) 
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Let us consider |p c ^), conventional proton wavefunction, as an example. Its 
spin-flavor wavefunction is (S is the symmetrization operator) 



b c t>SF = AfS[\D) nd <g> |«t>] = MS [( |ut di) - \ui dt) ® |«t)" 

= -AA £ »«A 5 H> n+1 Ml)™ H> nd ~ n Mt) nd " n , (35) 



where N c = 2nd + 1 and M is the normalization factor. For the £77(4) bases, 
it is written as a superposition of states with different R33 eigenvalues. Using 
the i? 33 eigenvalues for (u"\,d\.,u\.,d J {) = (1, 1,-1,-1), we can write 



n d 



|p c t)sF = JV'£ n d C n (-l) n I 1 , 1 , 4n + 1 - 2n d ) MS • (36) 

n=0 

This spin-flavor wavefunction contains states with both positive and negative 
R33 eigenvalues. The mean value and fluctuations of the R33 charge can be 
comparable during the time evolution. In addition, this makes the evaluation of 
i?33 expectation value for the conventional nucleon state much more nontrivial 
than for the hedgehog state. In fact, the contributions from states with positive 
and negative i?33 eigenvalues might strongly cancel one another. Our proposal 
of the dichotomous nucleon wavefunction, which will be discussed next, was 
motivated by this observation. 

The axial charge of the conventional nucleon wavefunction has been computed 
in seminal works, and is given by 

rid AT 1 9 

g\= (p c t^33|p c t)sF= |A/T ]T UA) 2 (4n + l-2n d ) = , (37) 

n=0 6 

where the contribution of 0(N C ) survives even after cancellations. Flipping 
either spin or isospin changes the sign of the axial charge, so g c A can be ±(JV C + 
2)/3 for nucleons. 



3.3 A dichotomous nucleon wavefunction 



The dichotomous wavefunction is constructed in such a way to make a nucleon 
axial charge g A small. To see how to achieve this, let us first consider the 
maximally symmetrized spin-flavor wavefunction made of = (N c — l)/2 
diquarks with I — S — 0, 

rid 

\D) nd ] = E n d C n (-1)" I , , An - 2n d ) MS • (38) 

n=0 
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If we take the expectation value of R33 for this state, we have 



-R33 <S 



\D) 



£ ( nd C n y (4n-2n d ) = 0, 



(39) 



n=0 



where the contributions to the axial charge completely cancel out. From this, 
we can see that in the conventional nucleon, the 0(N C ) contribution arises 
from the overlap between the diquark part and an unpaired quark. 

Therefore, if we place the unpaired quark in a spatial orbit different from other 
(N c — l)/2 quarks, the 0(N C ) contribution to the axial charge disappears. The 
dichotomous wavefunction for a proton with spin up is written as 



p d |) = MS \D;A) nd <g> \uf,B) 



(40) 



where A and B characterize spatial wavefunctions of quarks participating in 
the diquarks and an unpaired quark, respectively. When A and B are taken to 
be completely orthogonal, the axial charge is saturated by that of the unpaired 
quark, 



g d A = (p d t I R33 I /t) = (ut I ^33 I «t> = 1 • 



(41) 



So for the dichotomous nucleons, g\ = ±1. The more general case where A 
and B are not completely orthogonal was computed in Ref. [1]. 

The main problem when we try to regard this state as the ground state is the 
kinetic energy cost when placing an unpaired quark in an orbit different from 
the others. Its energetic cost is O(Aqcd)- 



4 Perturbative analyses 



Now we have prepared all ingredients necessary for our perturbative consid- 
erations. In this section, we compare the nucleon mass for different wavefunc- 
tions, at the level of the perturbative contributions up to 0(R 3 ) operators. 
It turns out that the hedgehog, conventional, and dichotomous wavefunctions 
give the same self-energy at the leading N c contribution. At NLO of N c , the 
signs flip for different orders of the i?-operators, so we can not conclude which 
wavefunctions give the lowest energy at this level of analyses. 

We will also argue specialities of the hedgehog state in these perturbative 
estimates, which will justifies the static mean field picture of the pion fields 
for the hedgehog quark wavefunction. 
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4-1 The contributions from the R 2 operator 



To compute the one pion exchange contributions proportional to R 2 ai , we use 
the relation 

Rl = C a - rl - a 2 . (42) 

where C 2 is the second order Casimir. The largest contribution to R 2 ai comes 
from C%. For the maximally symmetric spin-flavor wavefunction of n quarks, 
C*2 is given by 

C 2 (MS(n)) = 3n(n + A), (43) 

(see Appendix. B. 3). The hedgehog and conventional nucleon wavefunctions 
belong to the maximally symmetric spin-flavor representation of N c quarks, 
and the expectation value of C 2 is given by 3N C (N C + 4). 

On the other hand, in the dichotmous wavefunction, (N c — 1) quarks and an 
unpaired quark have no overlap so that the expectation value of the Casimir 
is given by a sum of that of the (N c — 1) quarks and one unpaired quarkp"} 
Then C 2 is given by 3(N C - 1)(JV C + 3) + 15 = 3N 2 + QN C + 6, which differs 
from other two wavefunctions by 0(N C ). 

The expectation values of t 2 and a 2 are trivial for the conventional and di- 
chotomous wavefunctions with definite isospin and spin quantum numbers. 
The values are (t 2 ) = (of) = 3. 

On the other hand, for the hedgehog state as a mixture of various isospin 
and spin states, we need a little manupulations. The hedgehog state contains 
isospins and spins of 0(N C ), so naively one might expect 0(N 2 ) contributions 
for (t 2 ) and (erf). It turns out, however, that (r 2 ) = (of) = 3iV c due to 
cancellations of 0(N 2 ) contributions (See Appendix.C.l). 

Assembling all these results, the one pion exchange contribution is propor- 
tional to 

N c (hedgehog) 
- (R 2 a .) = SN 2 - 6 x \ 2N C -1 (conventional) (44) 

N c (dichotomous) 

after multiplying the overall sign for the potential energy. We do not find the 
difference at the leading order of N c . The conventional wavefunction acquires 
the largest attractive force at 0(R 2 ). Since the intrinsic quark- meson coupling 
is g^qq ~ N~ 1 / 2 , the energy difference appears to be ~ g 2 qq xiV c ~ Aqcd- 



To see this calculation explicitly, the use of Eq.(B.22) may be helpful. 



18 



Through Eq.(42), one can observe several interesting consequences for the 
excited states in the maximally symmetric representation. First, the single 
pion exchanges produce the same mass splitting rule for the higher isospin 
and spin multiplets. Second, for baryons with larger the isospins and/or spins, 
the self-energy contribution from the pion exchange is less important. 

Finally we have to mention about the size of the 1/N C corrections. The ex- 
pression (44) explicitly contains the subleading terms of 1/N C , and one can 
observe that the l/N c corrections are substantial for (R 2 ai ). We will return to 
this issue at Sec. 6. 



4-2 The contributions from the R 3 operator 



To compute the _R 3 -operators, we will use the relation, 

tijktabcRaiRbjRck = — C 3 — 8C 2 + 6R ai r a ai , (45) 

where C 3 is the third order Casimir. The relation directly follows by rewriting 
the expression of C 3 (See Appendix. B. 4). As before, the largest contribution 
comes from C 3 which appears to be 0(N 3 ). For the maximally symmetric 
spin-flavor wavefunction of n quarks, C 3 is 

C 3 (MS(n)) = 24n(n + 4)(n-2). (46) 

For the hedgehog and conventional wavef unctions, we substitute n = N c , and 
we get (C 3 > = 2AN C (N C + A)(N C - 2). 

For the dichotomous wavefunction, we have a sum of C 3 's for (N c — l)-quarks 
and an unpaired quark. We get 24(iV c - 1)(JV C + 3)(iV c - 3) - 120 = 24(iV c 3 - 
iVp — 9iV c + 4), which differs from other two wavefunctions by O(N^). 

We already know the expectation values of C 2 . The remaining Rra term yields 



(Rai7~aO~i 



—N 2 

c 

3(iV c + 2) 
9 



(hedgehog) 

(conventional) 

(dichotomous) 



(47) 



The computation is given in Appendix. C.l. In particular, for the conventional 
and dichotomous wavefunctions, we will see (R a iT a cr,i) = 9|(i? 33 )|. 



Assembling all these results, the contribution from the uj — 3n coupling is 
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proportional to 



f 7iV c 2 + 8iV c (hedgehog) 
-e a bceijk{RaiRbjRck) = + 6 x j 6iV c 2 + 5N C - Q (conventional) (48) 

[ 3iV 2 -iV c + 3 (dichotomous) 

after multiplying the overall sign for the potential energy. The hedgehog state 
has the largest repulsive force. Again we do not find the difference at the 
leading order of N c . Considering the overall factors of the Wess-Zumino vertex 
and the quark number density of 0(N C ), the energy difference appears to be 
~ N c gl qq /ft x N c x iV 2 ~ A QCD . 

Again we can explicitly see the size of 1/N C corrections. The 1/N C corrections 
are substantial as seen in (-R 2 ,;). 



4-3 The hedgehog wavefunction and mean field pions 



The hedgehog state has special properties in processes with perturbative pions. 
For linear i?-operators, it has large axial charge distributions of 0(N C ), 

(H\R 33 \H) = {H\R+-\H) = {H\R+.\H) = -N c , (49) 

or (R11) = (R22) = (R33) = —N c . All the other expectation values of the R 
operators vanish. 

Let us first see the i? 2 -operator. As we have seen, the sum of the axial charge 
operator is ~ 3iV 2 . It means that the hedgehog state well-saturates the con- 
tributions from the intermediate states: 

E(H\Rlm = EE\(H\R m \n)\ 2 > j:\(H\R ri \H)\ 2 = 3N* C , (50) 

a,i a,i n i=l 

where we used the hermiticity of R a i. 

The similar relation also holds for the i? 3 -operator. Its leading iV c contribution 
is —6N^, which can be obtained by saturating the intermediate state with the 
hedgehog state 

£ijk£abc(H\RaiRbjRck\H) = 6{Ru)h{R22)h{R33)h-\ = ~ 6A,, 3 H , (51) 

where • • • gives the remaining 0(A" 2 ) contributions. 

As we have just seen, for the hedgehog state, the distribution of the axial 
charge operators may be treated classically at large N c . Correspondingly, the 
mean field picture for pions is at work for the hedgehog state. 
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In contrast, usual baryon wavefunctions do not have this factorization prop- 
erty, and in fact, large contributions come from the transition amplitudes from 
one baryon to other baryon states. Thus quantum treatments are essential for 
these wavefunctions. In particular, for the dichotomous wavefunction, even if 
the expectation value of the linear operator R (or the value of Qa) is small, 
the contributions from the off-diagnoal matrix element are very large. 



4-4 Color magnetic interaction 

Here we shall consider one gluon exchange type processes and see how it affects 
the energy splitting among different wavefunctions. 

The exchange of the spatial gluons (color-magnetic interaction) provides the 
spin splitting potential, (A = N c g^, and r qq i is distance between q-th. and g'-th 
quarks) 



which has been regarded as a source of the N — A splitting in the constituent 
quark models^] Assuming all constituent quark masses are the same and 
computing the spatial part (ip(r qq '): relative spatial wavefunction), the matrix 
element can be written as 



The second term is common for differernt wavefunctions, so the difference 
comes from the first term. For conventional (or dichotomous) wavefunction 
with fixed spin of 5* = 1/2, (of) = 3, while for the hedgehog state, (erf) = 3N C . 
Since the interaction is repulsive, the color magnetic interaction lifts up the 
mass of the hedgehog state by ~ Aq CD , compared to the conventional (or 
dichotomous) one. This is a consequence that the hedgehog state contains 
higher spin states. 



6 There is some objection to this reasoning. Liu and Dong performed a lattice 
simulation with and without the Z-diagrams which represent the meson exchanges. 
After removing the Z diagrams, the N — A and tt — p are energetically degenerate, 
while the one gluon exchange should not be affected. For details, see Ref.[33]. 




(52) 




(53) 



21 



5 Non-perturbative analyses: Chiral Quark Soliton 



So far we have considered the perturbative contributions, provided that our 
wavefunctions have the same constituent quark bases with the same quark 
mass and quark-meson coupling. Then we found that different wavefunctions 
give the same baryon self-energy at the leading N c . However, once a large 
pion field is developed, the quark bases can be modified, and then the 0(N®) 
modification for each quark mass may generate 0(N C ) difference for a baryon 
mass. We shall take into account such a possibility within the framework of 
the chiral quark soliton model. 

As we discussed in the previous section, the static classical picture of the pion 
field should work for the hedgehog state, while it is difficult to apply this 
framework for the conventional and dichotomous wavefunctions. Therefore in 
this section, we consider the hedgehog wavefunction only, and discuss whether 
the generation of coherent pions is favored or not. 



5.1 A baryon with the pion mean field 



The chiral quark soliton model is a model which smoothly interpolates the 
constituent quark model and topological soliton model. We will work in the 
Euclidean space. We start with the following Lagrangian 

C = qD(U 5 )q, (D(U 5 ) = (tf - MU 5 ) , U 5 = exp [i 75 r a 7r a ] ) , (54) 

for which the regularized partition function is 

Z = J VUVqVqe i $ d4xC | = J VU exp [iV c Tr Ln (d(U 5 )/D(U°)) ] , (55) 

where U = exp(ir a 7r a ). Here we divide the partition function with general U5 
by that with the vacuum configuration, U® = 1. 

For the vacuum, the derivative expansion around U® = 1 generates the non- 
linear a model plus an infinite number of higher derivative terms, that are 
powers of dU$/M, and the Wess-Zumino term [34]. The generation of the 
coherent pions should cost energy. 

The single quark bases for a given pion background is obtained by solving the 
eigenvalue problem (r = ix°) 

lo D(U 5 )<f> n = (-d T - h(U 5 j) . ( h(U 5 ) = 7o ( hkd k + MU 5 )) (56) 

We will denote the single particle energy for general U 5 and U® as e Q , e°, 
respectively. These bases are used for the evaluation of the quark determinant. 
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To compute the baryon spectra, we must consider the baryonic correlator at 
large time separation. It is given by 



(J B (6,T)J B (0,0)) U5 ~ JvUexp[N c TiLn(D(U 5 )/D(U° 5 ))} f[S q (0,T;U 5 ) 

i=l 

~ JvUexp[-T (^eafC/s] + E V 4U 5 })} , (57) 
where E sea comes from the regularized Dirac determinant [25], 

iV c TrLn (D{U 5 )/D(U°)) = —TN C £ (e a - e°) = -TE sca (U 5 ) , (58) 

e a ,e° a <0 

and is 0(N C ). This part is considered to be responsible for the polarization 
of the media that generates a pion cloud. Actually even after subtracting the 
trivial vacuum contribution, the sea contribution still diverges and has to be 
regularized by the UV cutoff 7 ] The coherent pion field affects quarks in the 
Dira sea with large negative energies. 

On the other hand, the E va x is the valence quark contributions originating 
from the interpolating fields (the product of iV c -propagators), 

-Evai = e vai(^5) > ( e vai(^s) : the level above the Dirac sea) (59) 

a=l 

which is also 0(N C ). Here the e va i must appear above the Dirac sea because 
of the Pauli-blocking. 

While the formal expression of the baryon mass is rather simple, its computa- 
tion is quite tough. The derivative expansion around the vacuum configuration 
is hard to control, because the pion configuration varies at the scale of Aqq D . 
It means that all the higher powers of dll^/M equally contribute, so that it 
is necessary to completely diagonalize the quark determinant as well as the 
valence quark propagator. Such treatments have been carried out numerically 
in seminal works [25,26,27]. We will borrow their results for our discussions. 



5.2 The optimized configuration within the stationary phase approximation 



Note that the valence quark contributions are 0(N C ), thus can affect the opti- 
mal configuration of U$ at large N c . The 11$ may deviate from 1 at the classical 
level. We determine the optimal configuration {7 5 from the following equation, 



In a more realistic treatment, there should be a form factor effect for the quark- 
pion coupling, which is likely to remove the sensitivity to the UV cutoff. 
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5u 5 (E sea [U 5 ]+E val [U 5 \) 



u 5 =ug 



0. 



(60) 



For two flavor case, we may parameterize U 5 as 



U 5 = exp ij 5 Q(x)n a (x)T t 



a 



cos 0(f) + vy^n a (x)T a sin 0(f) , 



(61) 



where = 



1. It is related to the linear realization as 



a(x) = M cos 0(f) , TT a (x) = n a (x)M sin 0(f) . 



(62) 



The 0(f) is responsible for the fraction of the chiral density and pion density, 
while n a (x) characterizes orientations of pions in the isospin space. 

In Eq.(60), if we take the variation with respect to 0(f), it yields 



Since the expectation value of the operators depend on the quark bases, this 
is the self-consistent equation. 

5. 3 The behaviors of E Rc i d and E va _\ 

Consider first the trivial configuration U5 — 1 (0 = 0) for which E sea = by 
definition. Then the baryon mass is determined by the lowest quark energy 
orbit without a pion cloud, and M B = N C M at this level of the treatment. 
Thus if a nontrivial configuration of U5 would exist, at least it must give 
smaller baryon mass than N C M. 

Now suppose that the hedgehog quark wavefunction would give the ground 
state. Through the Dirac equation for the quark hedgehog wavefunctions [36], 
the form of n a (x) is fixed to n a (x) = r a /r, and 0(f) = 0(r) where r = f — f 
(fo: the center of the baryon). 

For this configuration, the boundary condition at r — > is imposed to avoid 
the singularity related to the angular variable n a or n a . It requires 0(r) — > nn 
(n: integer) or ir a (r) — > as r — > 0. The integer n appears to correspond to 
the topological number of the meson profile. On the other hand 0(r) — > as 
r — > 00 to recover the vacuum configuration at asymptotic distance. 

We have to solve the self-consistent equation (63) numerically. Instead, we use 
the profile function 



sin 0(f) (qq(x) ) - n a (x) cos 0(f) (qi^r a q{x ) ) = . 



(63) 



Q(r;R) 



7rexp(— r/R) , 



(64) 
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R 



Fig. 6. A schematic plot of E VSj \, E SC3i , and their sum E to t as a function of a meson 
cloud size R in the profile (64). The decreasing behavior of E VSL \ is due to a larger 
binding energy for a broader potential well. The increasing behavior of E sea comes 
from the Dirac sea polarization. In -Etot) the minimum is typically found at some 
value of R between M _1 and 2M _1 , but the details about the height and locations 
of the minimum depend on the UV cutoff for the model. Here we took the UV cutoff 
to make the sea contribution small enough. With too large UV cutoff, the curve for 
E sea is lifted up, making the minimum a local one. 



to see the general tendency of a single particle spectrum and the Dirac sea 
contributions in the presence of the hedgehog background. Actually, the self- 
consistent solution of Eq. (63) is known to give the similar form [26]. Here R 
is treated variational parameter. 

Clearly E sea is a monotonously increasing function of M or R. Meson profiles 
with larger R and M just increases the energy cost to make a hole in the 
vacuum condensate. On the other hand, the single particle energy E va \ is a 
monotonously decreasing function of M or R. The couplings between quark 
and mesons are attractive, so larger M gives larger energy reduction. For larger 
R, the quark can be more deeply bound to the mesons, because of the smaller 
kinetic energy in the potential with broader range. At MR 3> 1, the valence 
energy level even enters the negative continuum, approaching to —M. 

These tendencies are summarized in Fig. 6 for our trial meson profile (64). 
The minimum (or local minimum) typically appear around R ~ (1 — 2)M _1 
after cancelling out very large binding energy and sea quark contributions. 
The corresponding baryon is a mixture of valence quarks and coherent pions. 

Actually whether the minimum appears below or above depends on subtle 
details in calculations, such as the UV cutoff. The larger UV cutoff increases 
E 8ea because in this model coherent pions even affect the energy level very 
deep in the Dirac sea, costing energy [26]. But in what follows, the existing 
literatures suggest that even if we find a nontrivial minimum with a coherent 
pion cloud, the energy reduction from N C M is 20 — 30% at most. 
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M 1 2M 1 



Fig. 7. An expectation when the confining effects are included for the previous plot 
in Fig. 6. (The confining potential lifts up the quark kinetic energy so even at R ~ 0, 
so the energy is higher than N C M.) The main modification appears in the behavior 
of Eygi- At distance scale r ~ M _1 , the linear potential closes a potential well made 
of a pion cloud, so the valence level does not acquire energy reduction much. As a 
consequence, the local minimum in E to t would disappear. 



5.4 Spatial size of coherent pion and confining effects 



The reduction of the baryon energy only comes from the quark binding to a 
pion cloud. The spatial size of the pion cloud must be so large, R > M -1 , that 
bound quark wavefunctions can spread out to avoid the kinetic energy cost. 

However, the quark wavefunction can not be widely spread out once confining 
effects are turned on. Our expectation is summarized in Fig. 7. The confining 
force is activated when a quark is at distance of r ~ Aq^d from the center of 
the baryon, giving the energetic cost of ~ ar ~ A QCD . When we add contri- 
butions from the pion cloud potential and the linear potential, the potential 
well in the former with depth ~ M is always closed at distance r ~ M' 1 by 
the linear potential, no matter how R is large. Then the deeply bound valence 
quarks may not be allowed even at very large R. 

As a consequence, it is possible that the nontrivial minimum, that was found 
within the chiral quark soliton model, disappears by including confining effects. 
As we have mentioned, the energy reduction by the quark binding to pions 
is not so large, so it is not surprising that confining effects lift the energy of 
the baryon with coherent pions above those without coherent pions. If this 
is indeed the case, the true energy minimum should be found at the small 
size of the coherent pion, R <C M~ x . Then a number of pions becomes small, 
invalidating the mean field or coherent pion picture. We have to go back to 
the quantum calculations to treat quantum pions rather than coherent pions. 
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6 Summary 



In this paper, we have discussed the energy differences among the hedgehog, 
conventional, and the dichotomous nucleon wavef unctions. Mainly we focus on 
the self-energy difference originating from the pion exchange which is related 
to the matrix element of the axial charge operators. 

After examining both the perturtative and non-perturbative treatments of 
pions, we realize that it is very hard to derive decisive conclusions on the lowest 
energy state. The perturbative treatments give the same leading iV c behaviors 
for our three wavefunctions. In the non-perturbative treatments within the 
chiral quark soliton models, the existence of the state with coherent pions 
seems to be a subtle issue and depends upon details such as the UV cutoff 
and the confining force which are difficult to include in a solid way. 

On the other hand, there seems to be a general tendency that both one gluon 
exchange and confining interaction generate extra energy costs for the hedge- 
hog wavefunction, which is a representative state with coherent pions. Both 
effects are not manifest in typical chiral effective models. The combination of 
these effects with chiral models deserves further studies. 

Clearly our analyses are incomplete. We dropped off exchanges of higher towers 
of mesons that couple to the axial charges. We expect that such resonances 
just change an overall size of the interactions, but do not alter the ordering of 
three wavefunctions. We have to check more carefully whether this expectation 
is correct or not, though. 

Another oversimplification is the nonrelativistic treatments. In particular, it 
is often said that at N c = 3, the nonrelativistic value of the axial charge is 5/3 
for the conventional nucleon, but it can be reduced to the experimental value 
~ 1.26 by including the down component of the Dirac spinor. 

However, presumably this is not the end of the whole story. The nucleon is 
measured on the lattice with several pion masses [37,38,39]. The value is typ- 
ically lower than the experimental value by about 10%. What is nontrivial is 
that even at a large pion mass ~ 700 MeV, g& remains small, ~ 1.1 — 1.2. The 
nonrelativistic treatments should work better for larger m n and higher con- 
stituent quark mass, reducing corrections from pion dynamics and relativistic 
effects. Thus we need explanations other than the relativistic corrections. 

As for the lattice studies, it should be also mentioned that a lattice study 
of the baryon spectra is recently carried out for iV c = 3, 5, 7, in which the 
lowest pion mass is (m^/m^) 2 ~ 0.2 [40]. The results show the rigid rotor 
spectra for the excited baryons, that are consistent with predictions from the 
color magnetic interaction, the perturbative pion exchange, and the collective 
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excitations of the chiral soliton. Again we can not judge which predictions 
are prefered, but we hope that eventually detailed studies of the pion mass 
dependence will distinguish these three contributions. 

Finally we have to admit that the l/N c corrections are substantial, at least 
in the perturbative estimates of the self-energy contributions (See, the ex- 
pressions, Eqs.(44) and (48)). This fact may suggest us to revise our original 
proposal (see, introduction), by taking a more appropriate leading order of 
the expansion. 

Maybe a part of the problem comes from the fact that we put a large number 
of quarks into a small flavor space (two-flavor). It is interesting to check what 
happens if we take the Veneziano limit, N c ,Nf — > oo, Nf/N c = fixed. For 
instance, the second Casimir for the maximally symmetric representation is 
given by 



It is clear that the the 1/N C correction appears as the Nf/N c correction, and 
is substantial. Other quantities and their qualitative impacts on the nuclear 
physics will be discussed elsewhere. 
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A More on the SU(4) algebra 

For later calculations, we wish to write commutation relations in terms of the 
Cartan bases [42]. We choose the following Cartan generators (Hi, H 2 , H$) = 
(r 3 , 03 , -R33). To find other bases, we must take appropriate linear combina- 
tions of r, a, and R to find generators E± a such that 



C 2 = (2N l -l)xN c (N c + 2N l ). 



(65) 



[ Hi , E± a ] — ± aiE± a . 



(A.l) 
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Let us find such linear combinations. First we define 



( 9 ) _ A q) ± ir 2 (g) (q) _ a? ± u#> 
T± ~ 2 ' a± ~ 2 ' 

d(<?) _ ,(«).(«) d(?) _ p(9) _ p(?) _ -(9)^(9) /A o\ 

Not all of these generators may be directly used as E± a . In particular, o~± and 
t± operators can not be used as raising and lowering operators for the irre- 
ducible representations of SU(A), because states with definite i? 33 eigenvalues 
are mixtures of states with different spins and isospins. It turns out that our 
bases are 

E± a (i) = R± 3 + r± , E± a (2) = R±3 — t± , E± a ( 3 ) = R 3 ± + o± , 

E ±ai 4) = R 3 ± - o± , E ±a ( 5 ) = 2R±± , E± a ( 6 ) = 2R± T . (A.3) 

Here the root vectors are 



a 



(i) 



(2,0,2), a< 2 > = (2,0,-2), a< 3 > = (0,2,2), 



5 (4) = (0, 2, -2) , 3< 5 > = (2,2,0), ^ = (2,-2,0), (A.4) 
which give the relation 

[E aUh E_ aU )] = a^H i . (A.5) 



B The evaluation of Casimir operators 

Here we briefly summarize the construction of the Casimir operators for the Lie 
algebra. To derive the general form, we follow the discussions of Okubo [41]. 
Then we use it to explicitly derive the formula for the maximally symmetric 
(MS) representations. 

B. 1 The general form 

We consider generators of the Lie algebra 

[X a ,X b } = f ab c X c . (B.l) 

The general form of the Casimir operator is given by (C p can differ from C p 
by an overall factor which will be fixed later) 

C p = g a ^- a " X ai X a2 ---X ap , (B.2) 
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where X ai , ■ ■ ■ X ap , are generators of the group, and 

Here X{ Xo } means a generator in a representation {A } which we shall call a 
reference representation. For {A }, we use the fundamental representation of 
SU(2N { ) with trf un d. [1] = 4, and we raise and lower the indices using g ab = 
2N{5 a b and g ab = (g~ v ) a b = ^ab/2N{. We sometimes drop off the subscript {Ao} 
in g ai '" a p to simplify equations. 

We wish to evaluate the value of the Casimir for any representations. Actually, 
once we compute the value of the quadratic Casimir for given representations, 
we can express the higher order Casimir using CVs. Below we briefly explain 
a main idea of the derivation. 

Let us consider a product space {A} ® {A }, and its decomposition into irre- 
ducible representations, 

N 

{A}®{A } = £©{A,}. (B.4) 

3=1 

In the RHS, there can be degeneracies of representations, for instance, it is 
possible that {A&} = {A;} for k ^ I. We wish to compute the Casimir for the 
representation {A}, using the fundamental representation for {A }. 

Consider a generator in the product space and its decomposition to irreducible 
representations, 

N 

X {A®A„} = xi x} (g) E {Xo} + E {x} ® X| Ao} = ]T X^Pj , (B.5) 

where Pj is a projection operator picking out the j-th representation space. 
Clearly J2j Pj = £ {A} ® ^ {A(,} . 

Now let us define the following operator 

Q\,\q = x { a x} ® ^{a o} = g ab x { a x} ® ^ {A()} , (B.6) 

which commutes with with all the generators in {A <8> Ao}, 

Q x , Xo , = . (B.7) 

Using this property, we can further show that the p-th order product 

(Qx-mY = X { a X} X£ } ■ ■ ■ (S> Xll o] XH Q} ■ ■ ■ X a { l o} , (B.8) 
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commutes with all the generators, and is the Casimir operator in the product 
space {A ® Ao}. Now taking the parital trace for {Ao}, we obtain g ai '" a p. Then 
the operator tr Ao [Q p ] gives the Casimir C P (X) for the representation {A}. 

The relation between the C P (X) and CVs follows from the property of the 
operator Q\-\ . We can rewrite Q\-x as a sum of C2 for the representations 
{Xj}, {A}, and {A }, 

Q x . M = g ab xW®X { b Xo} 



/ 



i; [{X a X b ) {mXo} - {X a X b ) {x} <g> - (X a X b ) {Xo} 



1 N 

.E[W-^)-W]P r (B.9) 



Below we will write £a ; a,a — ^(Aj) — C 2 (A) — C 2 (X ) /2. Then the p-th 
order product is 



AT 



(QX-M) P = ■ ■ ■ X^ ® X- } • • • X&y = E (^;A,A ) Pj ■ (B.10) 



Now a full trace of the operator {Q\-\ ) p for {A} and {Ao} gives the p-th order 
Casimir operator for {A} times the dimension of {A} (or the Casimir operator 
for {Ao} times the dimension of {Ao}). The full trace in RHS gives a sum of 
the dimension for each representaion {Xj}. Therefore we finally find that 

CM - E (6,;A,A ) P , (B.H) 



where d(X) is the dimension of the representation {A}. To compute C P (X), we 
need to know CVs of all the irreducible representations out of the product of 
{A} ® {A }. 



B.2 The Casimirs for the maximally symmetric representations 



In the following we will consider the Casimir for the MS representation for n- 
quarks. We need to compute the product of the MS(n) representation and the 
fundamental representation which gives only two irreducible representations. 
In the Young tableaux, 

■ • -rm <g> □ = ■ ■ - rrm © y 1 1 1 i - • -rm (B.12) 

In the RHS, the first one is the MS wavefunction of (n + l)-quarks, and the 
other is the mixed symmetric representation, which we will denote "mS" . 
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The dimension of the MS representation for n-quarks is given by 

By comparing both sides of the Young-tableaux, we can find the dimension of 
the mS wavefunction with (n + l)-quarks as 



rf(mS(n + 1)) = d(MS(n)) x 2N f - d(MS(n + 1)) 
n(n + 2Nf-l)\ 
~ (n + l)!(2iV f -2)! ' 

Since 

d(MS(n + l)) _ n + 2N { rf(mS(n+l)) _ n (2N i -l) 
d(MS(n)) n + 1 ' d(MS(n)) " + 1 



(B.14) 



(B.15) 



we find 



C p (MS(n)) = fa^)' + (6^,,)' • (B.16) 

The remaining task is the evaluation of the quadratic Casimir for the MS and 
mS representations. 



B.3 The second order Casimir for MS and mS representations 

The quadratic Casimir can be written in terms of the Cartan bases as follows: 

C 2 = E^ + £{£ Q 0) ,E_ aU) }. (B.17) 

i=l j=l 

This is related to C 2 = C 2 /4. (C 2 = X a X a = g ab X a X b = 5 ab X a X b /4.) We will 
evaluate it for the MS and mS representations with n-quarks. 

For the MS representation, we need to consider only the state with the highest 
weight, (\u t)) ra = \n,n,n). In the second sum, E ±a (j)\n,n,n) = for j = 
2, 4, 6, so only terms with j = 1,3,5 can contribute for this state. Since = 
(2,0,2), a { V = (0,2,2), and = (2,2,0), we have 

C 2 |n, n, n) = (Yin? + 4(r 3 + <t 3 + ^33)^ \n,n,n) = 3n(n + 4)\n,n,n) , 

(B.18) 
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where we rewrite {E a ,E_ a } = [E a ,E- a ] + 2E_ a E a , and use the relation 

[E aU) ,E_ aU) ] = a? ) H i . 

Next we consider the mS wavefunction with (n + l)-quarks. First we prepare 
the highest weight wavefunction of the mS representation which is orthogonal 
to the MS wavefunction with the same (T3, 03, -R33) values. Consider the MS 
wavefunction with the weight of the highest value minus one root, 

|n-l,ra-l,n + l) M s oc S [ K) n |<4)] , (B.19) 
which can be rewritten as 

S[\ut) n \di)] =S[\ut) n ] <g> \di) +nS [ |ut) n_ Vi>] ® K> 
= n! ( |n, n, n) MS ® |d|) + |n - 2, n - 2, n) M s ® |«t)) • (B.20) 

Then the mS wavefunction can be constructed as follows 

1 71 - 1,71 - 1,71+ l) mS 

= , (v^ ™, ™)ms <8> Ml) - \n - 2, 7i - 2, n) M s <8> |«T)) • ( B -21) 
V " + 1 

To evaluate the Casimir value, we decompose the generators as follows: 

(n \ ^ / n \ n 

EX (,) +X (,0 = £j#> +2£*W ) + (X(«')) 2 , (B.22) 
9=1 / \<?=1 / 9=1 

where the -X"^ 9 ') acts on the (n + l)-th quark. The square operators give the 
Casimir values for the MS-representation for n-quarks and a single quark. 
Nontrivial contributions come from the product of n-quark operators and a 
single quark operator, which can be rewritten as 

n n I 6 \ 

E = E HM* + E {E^% + E%E%1) . (B.23) 

9=1 g=i V i =1 / 

When these product part act on the state |n, n, u)ms <8> \dl), only ifj and £" a (5) 
parts give the nonzero contributions. We find 

n 

J2xi q) Xi q,) \n,n,n) MS ®\dl) 

9=1 

= -n \n, n, n) MS ® \d\) + 4^ |ra - 2, 71 - 2, n) MS ® |ut) , (B.24) 

where we use the relation \n, n, n) = 2y/n \n— 2, n— 2, n) and -E^) |g?4- 

) = 2 |wt)- Similar calculations are applied to \n — 2,n — 2,n) MS ® |^T)> f° r 
which only Hi and E a (5) survive. After summing these two contributions, one 
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finds 



J2X^X^\n-l,n-l,n + l) mS = -(n + A)\n-l,n-l,n + l) mS , (B.25) 



9=1 



\n - 1,71 - 1,71 + l)mS 
1 71 - 1,71 - 1,71 + l)mS 



so that C 2 (mS(n + 1)) is given by 

,9=1 / 9=1 

= [c 2 (MS(n)) -2(n + 4) + C 2 (MS(l))" 
= (3n + 7) (n + l) |n - 1, n - 1, n + l) mS . (B.26) 

In summary, the quardratic Casimir for the MS and mS representations are 

C 2 (MS(n + l)) =3(n + 5)(n + l) = 4C 2 (MS(n + 1)) , 
C 2 (mS(n + 1)) = (3n + 7)(n + 1) = 4(7 2 (mS(n + 1)) , 

so we have 



(B.27) 



£MS(n+i) = y [C 2 (MS(n + l)) -C 2 (MS(n)) - C 2 (l) 



£mS(n+i) = y [C 2 {mS(n + 1)) - C 2 (MS(n)) - C 2 (l) 



3n 

n + A 
4 



(B.28) 



Now we are ready to evaluate all the Casimirs for the MS(n) representation. 



B.4 The third order Casimir and its explicit form 



Using ^MS(n+i) and £ m s(n+i)j the cubic Casimir for the MS(ra) is 



C 3 (MS(n)) = 4 3 C 3 (MS(nj) = 4 

= 24n(n + 4)(n-2). 



/ \3 n + 4 / n3 3n 

(?MS(n+l) ) + (U(n+1) J 



n + 1 
(B.29) 



where we used C 3 = X a X b X c tr[X a X b X c ] = X a X b X c ti[X a X b X c }/ '4 3 = C 3 /4 3 . 

The contents of the cubic Casimir operator are the following. Apparently, 
there are R 3 , r 3 , a 3 , R 2 r, R 2 a, Rr 2 , Ra 2 , Rra, r 2 a, to 2 terms, but nonzero 
traces come only from R 3 , r 3 , a 3 , R 2 t, R 2 a, and Rra terms. We write 

C 3 = C{R 3 ) + C(r 3 ) + C{a 3 ) + C(R 2 t) + C(R 2 a) + C(Rra) . (B.30) 

C(R 3 ) is given by 



C(R 3 ) — tr [R ai R b jR c k\ R a iRbjRck — — 4 ^abc^ijkRaiRbjRck ■ 



(B.31) 
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C(r 3 ) 

C(r 3 ) = tr [T a T b T c ] T a r b T c = 4i e abc r a r b r c = -8 r 2 , (B.32) 
and C(i? 2 r) is 

C(R 2 r) = tr [i? a ji? bj r c ] R ai R b jT c + tr [i? ai T c i2&j] R ai T c R b j + tr [T c i£ a i-R&j] r c R ai R b j 
= 4ie abc (^R ai R b jT c — R ai T c Rbj + T c R ai R h ^j 
= 4:ie abc R ai [R b j, t c ) + 2ie abc T c [R ai , R b j] 

= -8 (2^. + 3 r a 2 ) . (B.33) 
We can obtain C(a 3 ) and C(R 2 a) terms as well. Finally C(Rra) term is 

C(Rra) = tr [R ai T b <jj} RaiUaj + tr [RaiVjn] RaiCjU + tr \T b R ai a 3 \ T b R ai aj 
+ tr [ajRaiTh] <JjR ai T b + tr [r^-i^] T b ajR ai + tr [cr,-r 6 ii! a j] (TjT b R ai 

= 4 (2R ai T a <Ji + iTaOiRai + T a R ai Oi + OiR ai T^ 

= 24R at T a a t , (B.34) 

where we used [i? a j, r a ] = [i? ai , <Xj] = in the last line to change the ordering. 
Assembling all these terms, we have 

C3 = —Ae abc eijkR a iRbjRck — 32C*2 + 24R ai T a ai . (B.35) 



C Computations of various matrix elements 

C. 1 r 2 and o\ for the hedgehog state 



We compute r 2 and of for the hedgehog wavefucntion. For the third compo- 
nent of the isospin, we have 

1 N c 

(H\r 2 3 \H) = — J2 N c C n (N c - 2nf = N c . (C.l) 

Z C n=0 

For the other isospin operators, we have 

(H\ (r 2 + r|) \H) = l - (H\ {E aW - E aW , E_ aW - E_ a{2) } \H) 

= i (H\ (E_ aW E aW + E ai2) E_ a(2) ) \H) , (C.2) 

8 Note that the trace runs over SU{A) fundamental representations so a factor 2 
should be multiplied to the results of the SU{2) fundamental representation. 
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where we used the fact that E_ a (i)\H) = E a (2)\H) = (H\E_ a (2) = (H\E a (i) = 
0. (This is because = (2,0,2) = (2,0,-2), and the hedgehog state 
has the minimum of the axial charge, -R33 = —N c .) Let us further note that 

(H\ (E_ ai i)E a (i) + E a (2)E_ ai 2)) \H) 
= (H\ ([E_ aW ,E aW ] + [E ai 2 h E_ a( 2 } }) \H) = (H\(cP> - d^H^H) 
= -A{H\R 33 \H) =AN C . (C.3) 

Applying the same manupulation to the spin operators, the hedgehog state 
has the following expectation values, 

(H\r 2 a \H) = 3N c , (H\a*\H) = 3N c . (C.4) 
C.2 The expectation value of the Rra operator 

Next we compute R ai r a ai in C 3 for various wavef unctions. First we write it as 

RaiT a O-i = R 33 T 3 (T 3 + 2 (R +3 T_(T 3 + R^ 3 T + a 3 + R 3+ T 3 (T- + R 3 -T 3 CT + ) 

+ 4 (R ++ T-a- + J R_+r+o-_ + i? + _r_o- + + j R__t + ct + ) . (C.5) 

We first consider the conventional and dichotomous wavefunctions for the state 
\pt)- The isospin and spin parts can be computed readily, and we have 

(ptlRaiTaVM) =9(pt|#33|pt}, (C.6) 

where we used (p-[\R++\ni) = (pt\R+ 3 \nt) = (ptl-Rs+b^) = (ptl-Rrabt), 
and r_|pt) = bt)> etc. 

For the hedgehog wavefunction, it is useful to note that (r^ + a\ q ^)\H) = 
for each q. Then, for instance, 

(H\R 33 r 3 a 3 \H) = -N c (H\r 3 a 3 \H) = N c (H\r? o \H) = N 2 C , (C.7) 

or another example is 

(H\R ++ T-<J-\H) = -J2(H\(r^) 2 r^\H)=0, (C.8) 

q 

since is the fundamental representation. Doing similar manupulations, 
finally we arrive at 

(HlRaiTadlH) = -iV c 2 . (C.9) 
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